Interaction induced edge channel equilibration 



Anders Mathias LundcjB Simon E. Nigg, and Markus Biittiker 
Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve 4, Switzerland 

(Dated: October 14, 2009) 

The electronic distribution functions of two Coulomb coupled chiral edge states forming a quasi- 
ID system with broken translation invariance are found using the equation of motion approach. We 
find that relaxation and thereby energy exchange between the two edge states is determined by the 
shot noise of the edge states generated at a quantum point contact (QPC). In close vicinity to the 
QPC, we derive analytic expressions for the distribution functions. We further give an iterative 
procedure with which we can compute numerically the distribution functions arbitrarily far away 
from the QPC. Our results are compared with recent experiments of Le Sueur et al.. 
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Two decades ago, edge states (ES's) [l| were demon- 
strated to be a physical reality by creating a non- 
equilibrium population Q through selective injection and 
detection of carriers in different states along the same 
edge Experiments revealed that the inter-edge 

carrier scattering could be strongly suppressed 0, 0, Q 
over distances of 80 /im. Now in a scries of novel experi- 
ments the group of Pierre @, [tI has investigated the non- 
equilibrium distribution function in an ES as it evolves 
along a channel away from a QPC at which it is initially 
created. The experiments are carried out in a high mobil- 
ity two-dimensional electron gas at a filling factor ly = 2 
such that there is an outer (spin up) non-equilibrium ES 
and an inner (spin down) equilibrium ES. The distribu- 
tion function is measured with the help of a quantum 
dot (QD) sufficiently small to provide transmission only 
through a single resonant level, see fig. [1] The QD serves 
as an energy spectrometer and permits the reconstruc- 
tion of the distribution function in the outer ES. 

The experiments reveal two surprising features: First, 
the initial non-equilibrium distribution created at the 
QPC and calculated form non-interacting scattering the- 
ory differs only weakly from the measured one over dis- 
tances of close to one micrometer 6, 7j. At large dis- 
tances from the QPC, due to the Coulomb interaction be- 
tween carriers in the two ES's, the distribution function 
evolves into an equilibrium distribution function at an ef- 
fective electrochemical potential and temperature. The 
outer non-equilibrium ES transfers part of its energy to 
the inner ES. The two ES's equilibrate towards the same 
equilibrium distribution with the same temperature (but 
still at different electrochemical potentials due to lack of 
particle exchange between the two ES's). The second 
surprise of the experiments is the fact that the temper- 
ature of the distribution functions at large distance in 
the two ES's is lower than dictated by equilibrium ther- 
modynamic arguments The first surprise shows that 
relaxation due to inter-ES interaction is weak. The sec- 
ond surprise implies that equilibration occurs not only 
between the inner and outer ES's but that there must be 
an additional equilibration mechanism which cools the 



two ES's below what would be expected from inter-ES 
coupling alone. We propose that additional excitations 
in the bulk Q , which couple predominantly to the inner 
ES, have to be considered to understand this effect. Al- 
though the nature of these excitations remains unclear, 
the experimental findings of 0] are consistent with this 
hypothesis. For example it is found that when the inner 
ES is forced to form a short closed loop, then relaxation 
in the outer ES is strongly suppressed. 
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FIG. 1: (color online), (a) The experimental setup to mea- 
sure the electronic distribution function of an ES. The full 
(red) curve represents the measured outer ES while the dashed 
(blue) curve represents a co-propagating inner ES. The two 
ES's exchange energy via Coulomb interaction between x = 
and X = L. The initial distribution functions (b) relax, via 
energy-conserving particle-hole excitation processes, toward 
Fermi functions (c). 

The physics of ES's is often discussed within the frame- 
work of bosonization theory, where the elementary exci- 
tations have bosonic character and are of collective na- 
ture In contrast, we take the weak equilibration 
seen at distances of less than a micrometer as the start- 
ing point of a discussion which treats inter-ES interaction 
perturbatively ll| . The interaction is described in terms 
of two-body collisions. We use the equation of motion ap- 
proach for second quantized operators to derive an evo- 
lution equation for the distribution functions which re- 
sembles a Boltzmann collision term with the added com- 
plication that there are two different initial distributions 
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(one for each ES). Alternatively the Coulomb matrix ele- 
ments which appear in this theory can be taken from an 
RPA theory [l^] in which the electron densities in each 
channel fluctuate and interact through an (effective) ca- 
pacitance. To treat equilibration at longer distances we 
iterate numerically the solution for short distances. At 
large distances the distribution functions approach their 
equilibrium form dictated by entropy maximization. 

We describe the ES's in terms of scattering states 
XaE^x) with energy E and label a = o,i {i : inner, o : 
outer). The inter-ES interaction is given by 

Hint = \Yl I dEdE'U^{E',E)alE,a^E, (1) 

a 

where aj^^ ('^qb) is the creation (annihilation) opera- 
tor for the scattering state XaE and Ua{E' , E) is the 
potential operator for scattering a particle from E to 
E' in the ES a at the expense of a particle scatter- 
ing in the opposite ES a. Explicitly UaiEy , Ei,t) = 
J dE2dE2'V^^E^, E^E^alE^^{t)ag^j^^{t) in the Heisenberg 

picture and V^^^a' £^1^2 inter-ES electron-electron 

interaction matrix element for the scattering process 
{aEi,aE2) {aEy , Q.E2') ■ Using the Heisenberg equa- 
tion of motion ihdta^^{t) = [a^^[t),H], the electronic 
distribution function fa{E) in ES a can be found by 
evaluating (a|,^(t)a„s, (t)) = S{E - E')f^{E). The 
non-interacting distributions are fi{E) = f^.{E) and 
f^{E) = nf^^iE) +Tf^^{E), where /" ^ {l+'exp[(£;- 
IJL)/k^T]] ^ and T (7?.) is the transmission (reflection) 
probability of the QPC, see Fig. [TJ The chemical po- 
tential of the inner ES fii can experimentally be tuned 
independently of /xi and ^2 by using an additional QPC 
(not shown in Fig. [T]). To second order in the interac- 
tion matrix element the distribution is fa^ ~ +6fa\ 
where (see Supplementary Material for details) 
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fjE + huj)[l - f^{E)]S^ij^suJE,E + huj,uj) 



/°(i?)[l-/^(i? + M]^, 



su^su. 



{E + hLu,E,uj) 



(2) 



The first term contains the absorption potential fluctu- 
ation spectrum S'^jj^^jj^{E,E' ,uj) describing an ab- 
sorption of energy hio by the ES a while the ES a goes 
from energy E' to E. Likewise the second term with 
the emission fluctuation spectrum S\ 



describes the 

emission of energy fiw from the ES a to the ES a, which 
consequently leads to the transition E ^ E + Tiw m. a. 
The fluctuation spectra are to lowest order in the inter- 
action and defined by 27r(5(ti' + uj')S'^jj {E',E,uj) = 
{SU^{E,E',Ljy^^SUa,{E\E,uj'y^^), where JC/U) ^ 



the deviation from the average potential to first order 
in the interaction. The emission spectrum is found by 
interchanging the two SUa in the absorption spectrum or 
equivalently by changing the sign of lu. Explicitly, the 
spectra are found to be 

h JdE"\V^f%„^^^^^^„\'fjE")[l - f^{E" + fuv)l (3a) 
SsUcwS^' iE,uj) ^ 

h ldE"\V^r^„^^^„_,^j'f^iE" + hu.yi - aE'% (3b) 

where the interpretation in terms of emission and absorp- 
tion spectra is clear. By inserting these into Eq. ([2]) the 
similarity with the collision integral in the Boltzmann 
equation becomes evident. 

(2) 

Next we wish to calculate 6 fa {E). To this end, the 
inter-ES scattering process {aEi,aE2) {aEi/ ,aE2') 
needs to be considered. If the ES's are perfectly trans- 
lation invariant, then energy and momentum conser- 
vation together reduce the available one dimensional 
phase space enormously compared to higher dimensions 
[14j . This leads us to consider the more realistic non- 
translation invariant case caused by the fact that the 
ES's follow the equipotential lines created by the sam- 
ple edges and the impurity potential. Including this non- 
translation invariant ES physics leads to the presence of 
non-momentum conserving scattering pro cesses increas- 
ing the phase space substantially [15|, |l6[. The broken 
translation invariance is included into the model of the 
inter-ES interaction matrix element \V§"^^^ BiBal^' 
avoid modeling a specific geometry we perform a statis- 
tical average over the geometry of the ES's and thereby 
introduce the momentum breaking correlation length £p, 
which is smaller than the size of the region of relax- 
ation L. For simplicity, an effective interaction of the 
form V{x,x') = S{x — x')g{x) is used and it is as- 
sumed that the deviation of g{x) from some mean value 
(7o is Gaussian distributed, i.e. {g{x) — go){g{x') — go) = 
A/{V2^ep) exp [~{x~ a;') V(2£2)] where A/{V2^£p) is 
the maximal deviation and denotes the geometrical 
averaging. This yields an interaction with a momentum 
conserving and a momentum breaking part. The latter 
is (see Supplementary Material for details) 
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■exp[-(AHp)V2] , (4) 



(1) 



(Ua^) is the Fourier transformed operator for 



where Ak = {Ei - Ey)/{nva) + {E2 - E2')/{hva), using 
linear dispersion relations with different velocities Va for 
the two ES's. Note that for linear dispersions with dif- 
ferent velocities there is no phase space for scattering in 
the momentum conserving limit, Afc = 0, but in the very 
special (almost pathological) case = Va, momentum 
and energy conservation are equivalent leading to plenty 
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of phase space. The specific model for the interaction and 
the matrix element is not of great importance as long as 
it includes the physics leading to non-momentum con- 
serving processes, which in turn introduces a new length 
scale 

For energy conserving scattering, the model matrix el- 
ement Eq. only depends on the transferred energy in 
the scattering [17|, since Afc = uj{l/va — 1/wa)- This 
means that the energy integral in the fluctuation spec- 
tra of Eqs. ([3]) can be done analytically upon which it 
becomes evident that 5f^\E) cx T(l - T). Thus the 
greater the shot noise of the QPC, the faster the re- 
laxation is. The elementary scattering processes lead- 
ing to relaxation consist of a particle loosing energy in 
the noisy outer ES and a particle gaining energy in the 
noiseless inner ES as illustrated on Fig. [U (b). The 
matrix element introduces a new energy scale /S.E = 
{h/ip)vaVa/{va — Va), which limits the possible amount 
of energy transferred between the two ES's in the scat- 
tering process since the matrix element is proportional 
to e-e^'^/'^-^)'. In the limit that fcgT, |^2 - < \^E\ 
the distribution functions for the inner and outer ES can 
be found analytically to be 



^/P(i?) = -7'T(i-r)(M2-m) 



E~-{fii+fi2) , (5) 



(6) 



X \ [fl (E) /°- (E)] [ink,Tf + {E- ^-f] 



where 7^ = AL/[h'^vl^vl] and /ij±(^2-/Ui) is 

the maximal and minimal energy of particles affected by 
the scattering process in the inner ES. Here it is seen that 
the maximal available energy (apart from thermal exci- 
tations of order fc^T) is given by the energy difference 
M2 ~ Ml creating the step distribution. The scattering 
processes create a linear slope on the plateau of the dis- 
tribution of the noisy outer ES as shown in Figl^l The 
slope crosses the middle of the plateau and it is propor- 
tional to the noise of the QPC and the energy available 
/i2 — Ml- The inner noiseless distribution gets a tail on 
both sides of the Fermi level, which extends over the 
length of the plateau M2 ~ Mi- the general case, the 
distribution functions can be found numerically and the 
matrix elements in Eq. ^ have to be included in the 
calculation, but the transferred energy is still limited by 
AE. 

The above perturbative results apply for a short dis- 
tance L after the QPC and express the distribution func- 
tions at L in terms of the (unperturbed) distribution 
functions at the origin. Once the distribution functions 
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FIG. 2: (color online). Analytically calculated inner (left) and 
outer (right) ES distribution functions in the regime k^T, \ fJ,2 — 
fill <C |Ai5|. The parameters are (energies in /xeV): T — 0.5, 
Ml = 15, M2 = 30, fii = 20, fceT = 0.2. 



at L are known we can use them to calculate the distri- 
bution functions at a distance 2L via Eq. (O. By iter- 
ating this procedure we can thus describe the effective 
length dependence of the energy relaxation. A conve- 
nient quantity with which to characterize the relaxation 
of fa at temperature T is given by the excess temperature 



defined as 



k T 

exc,a 



J dEAUE)iE - fia) - {k^Tf 



Here Afa{E) ~ fa{E) — 0{fla — E) is the difference be- 
tween the actual distribution function and a zero temper- 
ature Fermi distribution with the same number of parti- 
cles and hence /ia = £'0 + '^Efa{E), where Eq is cho- 
sen such that fa{E) = 1 for E < Eq. k-^Texc,a gives the 
energy of the non-thermal excitations in /„. The initial 
excess temperature right after the QPC of the inner ES 
is zero and the one of the outer ES is given by k^T^^^ ^ = 
{^T{1 — T)}^/^|/i2 — Mil- Because of energy conserva- 
tion, Tf,xc,a is a conserved quantity in the equilibra- 
tion process. Furthermore due to entropy maximization 
the excess energy is distributed equally among the two 
ES's, which in the limit of long distances thus converge 
towards Fermi distributions with equal excess tempera- 



tures given by k^T't 
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{^r(i - T)} 
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M2 - MiL- 



The excess temperature of the outer ES measured in 
does indeed saturate at large distances toward a finite 
value. This value is however found to be systemati- 
cally lower than the above prediction, for large volt- 
age biases. Surprisingly it agrees well with the value 



{T(l — T)}^/^|/i2 — Mil/"" expected from en- 



, ^00(3) 

"-B-*^ (^XC 

ergy equipartition among three instead of only two chan- 
nels. What could provide the additional relaxation chan- 
nel? In 01 it has been observed that if the inner ES is 
forced to form a short enough closed loop, such that the 
energy level spacing of its (discrete) spectrum is larger 
than the available energy provided by the voltage bias 
M2 ^ Mil then relaxation of the outer ES is strongly sup- 
pressed. Motivated by this observation, we suggest that 
there exist excitations in the bulk which arc coupled 
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via long range Coulomb interaction to both the inner ES 
and an ES on the opposite side of the sample. As long as 
the bulk excitations can be created, such a mechanism 
would allow extra energy to be carried away from the 
outer ES. 

As a first approach we model this extra degree of free- 
dom as an additional ES coupled to the inner ES only, ini- 
tially in equilibrium at the electronic temperature, which 
we take to be T = 30 mK. This then contributes an ex- 
tra collision term in Eq. ^ and allows a quantitative 
comparison with the experiment Q. The fitting proce- 
dure is detailed in Supplementary Material and the re- 
sult is shown in Fig. [H The best fit is obtained when 
the coupling strength to the bulk excitations is about 
three times larger than the inter-ES coupling strength 
and when AE = 14.3 /xeV, which for Vo and Vi between 
10^ and 10^ m/s leads to ip > 0.5 fiui. For intermediate 
distances (i.e. 2.2 //m and 4/zm) both the data and our 
numerics display a similar weakly non-linear behavior of 
the excess temperature as a function of the voltage bias. 



turn Hall regime. We derived an analytic expression for 
the leading order correction to the non-interacting the- 
ory, which is present in a system without translation in- 
variance. We have shown further that the result obtained 
in the long distance limit, by iterating the perturbative 
solution numerically, is in quantitative agreement with 
a recently performed experiment if we take into account 
the electric coupling between the innermost ES and bulk 
excitations. Finally we note that measurements of the 
distribution functions of both channels for even shorter 
distances than 0.8 /im should allow for further critical 
testing of our theory. 

The authors would like to thank H. Le Sueur, C. Al- 
timiras and F. Pierre for fruitful discussion and for shar- 
ing their data. This work was supported by the Swiss 
NSF. 

Note added. — During completion of this work we be- 
came aware of related work based on plasmon scattering 
by P. Degiovanni et al. (loj . 
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FIG. 3: (color online). The excess temperature of the outer 
ES versus the voltage difference across the QPC. For clar- 
ity curves for different voltages have been shifted upward 
by 75 mK. The dashed (black) lines indicate the initial 

value T^xc,o (upper lines) and the asymptotic value T^c'^^ 
(lower lines) as expected from energy equipartition. The 
thin (black) lines are guides for the eyes to better see the 
weak non-linearity present at intermediate distances (espe- 
cially for 2.2 /jm and 4/xm). Experimental data (circles with 
error bars), courtesy of F. Pierre et al.. 

To conclude, we have computed the distribution func- 
tions of two Coulomb-coupled ES's in the integer quan- 
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Supplementary Material for 
"Interaction induced edge channel equilibration" 

111 this Supplementary Material we present details of the derivation of the leading order correction to 
the distribution functions as given in Eq. (2) of the main text. We also describe the model used for the 
Coulomb matrix elements given in Eq. (4) of the main text. Finally, we explain the iteration procedure 
and how we fitted our theoretical predictions with the experiment of 0]. 



THE EQUATION OF MOTION APPROACH AND THE DISTRIBUTION FUNCTION 

Our starting point is the Hamiltonian describing the dynamics of the electrons in the inner and outer ES's after 
the QPC: 

H= J2 /'dxV'l(x)T"(x)7A„(x) + i^ /'dxd.T>t(^)^t(^/)y^^(^^^/)^^ (8) 

a—i,o af3 

The first term describes the kinetic energy plus the single particle potential and the second term the inter and intra 
ES Coulomb interaction of the outer (o) and inner (i) ES's. The intra-ES interaction typically leads only to a small 
contribution to the relaxation due to the presence of both the direct and the exchange term as we shall see shortly. 
Therefore the focus in the paper is on the inter-ES interaction Vaa, using the shorthand notation a = SiaO + 5oai (for 
the opposite ES of a). Introducing the scattering state representation 

OqE = J dx X*aEi^)^a{x) <=> ^Jaix) = J dE XaE{x)aa{E) , (9) 

makes the single-particle part of the Hamiltonian diagonal. In this representation, the equation of motion for the 
annihilation operators ihdtaaE{t) = [aasit), H] in the Hcisenberg picture (i.e. A{t) = e'^*Ae~*^*) becomes 

ih—a^E{t) = EaaE{t) + / dE'Ua{E,E',t)a^E'{t). (10) 



dt 

In the case of inter-ES interactions only, the potential operator is given by 

Ua.{E,E',t)= J dE2dE2'Vg%,^E'EAEjtWE,it), (11) 

where , , defined in Eq. below, is the Coulomb matrix element for a transition from an energy E' to E 
in ES a and a simultaneous transition from energy E2 to E2' in ES a. A standard perturbation treatment to second 
order in V^" leads to the result 



5{E-E')f^^\E)^{ai^a^E,t ' ^ 5{E - E') [^fl{E) + 5f^^\E)) , (12) 
with the inter-ES relaxation given by 

5fl^\E) = {27Tfn J dc.di?'|Fs"|,+,^,^+,^s,|2 (13) 

X 'fjE + h^)[l ,CiE)]aE')[l - fl{E' + r^)] - flm^ - fl{E + nu:)]fl{E' + fi^)[l - fl{E')]_ . 

The combination of Fermi functions which appears here ensures the Pauli exclusion principle. Furthermore one can 
easily show that 

{5U^{E,E',uj)^''^6U^{E\E,uT^) = {2^fh5{oj + J) f di?2H^s"'l,+r^.,iJijJ V°(i?2')[l " /^(i?2' + M] , (14) 



where SUa{E , E' , ajY^^ = Ua{E , E' , uj)'-^^ - {Ua{E, E' ,uj)''^'^) and Uc,{E , E' , uj)'-^'^ is the Fourier transform of the first 
order expansion of pi|l . This then leads immediately to Eq. (2) of the main text. 
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The relaxation due to intra-ES interactions can be found in the same way. The result is as in Eq. (|T3)) with the 
important replacement: 

\^EE'+huj,E+hujE'\'^ ^ O I ^EE' + huj.E+huE' ^ ^E'+hu)E ,E+hu)E' 1^7 (15) 

direct term exchange term 

where the two final states were exchanged in the second term called the exchange term. The appearance of both 
the direct and the exchange interaction is due to the fact that the electrons are in the same conductor and hence 
quantum mechanical exchange processes are important. For the simple matrix element model used below, the direct 
and exchange elements exactly cancel. Although this cancellation is model dependent, the intra-ES contribution is 
strongly suppressed as a general feature and is therefore safely neglected in the present work. 

Finally, we note that the generalization to more than two coupled ES's is straightforward and formally results in 
additional collision terms in Eq. (fTS]). In the following, we make use of this generalization to incorporate the effect of 
an extra relaxation channel. 



MODEL FOR THE COULOMB MATRIX ELEMENTS 



The interaction matrix clement is 

VE',,%,,E^E2 = J da;da;'x;^£;^,(a;)x*^£;^,(a;')VaiC2(a;,a:')Xa,B,(a;')Xai£;i(a;), (16) 

where Va^a^ (a;, x') is the effective Coulomb interaction between an electron at x in ES ai and an electron at x' in ES 
a2- In the adiabatic edge channel description jT3|, ES's follow the equipotential lines in the sample. Due to impurities 
near the edge and to the variation of the confining potential, the inter-ES distance is not constant in the direction 
of propagation. In other words, translation invariance is broken and momentum is not conserved in the collision 
process. For the sake of generality, we do not model a specific geometry nor impurity distribution, but instead simply 
average over the geometry. This naturally introduces a momentum breaking correlation length scale (.p quantifying 
the amount of momentum breaking present in the system. 



i 
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FIG. 4: (color online). Schematics of 
two co-propagating edge states with 
varying inter-edge state distance. 



As a specific simple example, we consider the interaction between points with the same x coordinate as illustrated in 
Fig- m i-G- Vaaix,x') = S(x — x')g{x) , where g{x) is an unknown function of the interaction strength variation along 
the ES's. This approximation is justified for very good screening. However, it is important to emphasize [l^ that 
this simplification is not of vital importance for the physics discussed here, i.e. if we consider the poor screening limit, 
then the same physical effect of the broken translation invariance remains as can be shown from a more complete 
calculation. We assume linear dispersion relations with constant velocities Va and use a single particle basis of plane 
wave states 

Xc.is(x) = ^^e»^°(^)^ with fc„(iJ) = A. (17) 
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Note that these states are normahzed to have equal current and obey / da;x*^(x)xQE' (a;) = S{E — E'). With the 
above assumptions, we obtain 

VS:,E^„E.E. = dxe^^^^g{x) , (18) 

where the amount of broken momentum conservation Afc = ka{Ei) + fcg(i?2) — kaiEy) — A;g(£'2') is introduced and 
X is integrated over the region of relaxation of length L. Next the geometrical averaging of the squared matrix 
element \V^"^^^ BiB, P is performed by assuming that the deviation Ag{x) = g{x) — go is Gaussian distributed, i.e. 

Ag{x)Ag{x') = A/{\^2tt£p) cxp[~{x — x')^/(2£p)], as described in the main text. This gives two contributions to 

\^e"e , E-iE2\'^'' with Afc = and one where Afc 7^ is possible. The latter is given by Eq. (4) of the main text. 
Note that it is proportional to L, because interaction is included only over a region of length L. 

In this work, linear dispersions with different velocities are used, which leaves no phase space for scattering in the 
momentum conserving case Afc = due to simultaneous energy and momentum conservation. In general, the disper- 
sions are not linear and this changes the phase space constraints due to momentum conservation in the translation 
invariant limit. For instance, for a spin-split quadratic dispersion there would only be a single possible final state, 
where the two initial momenta are simply interchanged. This would lead to a very small resonant-like feature on the 
plateau of the outer ES distribution function, which could be tuned along the plateau by the chemical potential of 
the inner ES /x,;. Any dispersion with a positive curvature will lead to the same result. However, this is a very small 
increase of the phase space compared to the increase introduced by the non-translation invariant physics and impor- 
tantly, this is not what is observed experimentally. Therefore for the non-translation invariant physics the curvature 
of the dispersion is without importance and for simplicity the dispersion is linearized. 



NUMERICAL ITERATION PROCEDURE 



Eqs. (2) and (3a-b) of the main text express the variation of the distribution functions of the ES's after a short 
distance L as a functional of the distribution functions at the origin. We thus have a system of coupled equations for 
the approximate distribution functions after a distance L, which is valid to second order in the interaction and has 
the form 

fo..L{E)^F[f^,oJo..»\ + 0{V^). (19) 
Iterating this, we obtain a recursive relation 

fa,nL{E) « F[fa^(n-1)L, /a,(n-l)i] > (20) 

where n is the iteration number. Thus wc can calculate an approximation to the distribution functions at distances nL. 

A typical result of this iterative procedure for three coupled ES's is shown in Fig.[5l where we show the evolution of 
the excess temperature of the inner ES, the outer ES and the additional relaxation channel (denoted with "bulk") as 
a function of the iteration number. For the numerical integration, we use a three points Simpson extrapolation rule. 
We see that the inner ES and the additional relaxation channel, which are initially at equilibrium both gain energy 
since their excess temperatures increase. This increase happens at the expense of the initially out of equilibrium 
outer ES, which loses energy during the relaxation processes and thereby its excess temperature decreases. All three 
converge towards the temperature ksT^c^^ = {T(l — T)}^/^|/ii — /i2|/7r, expected from energy equi-partition. Inset 
(a) shows the distribution functions after 1 iteration and inset (b) after 54 iterations. Note in particular that all three 
distribution functions are found to converge towards Fermi functions (dashed (black) curves in inset (b) which are 
indistinguishable from the calculated distributions). While the chemical potentials of the equilibrated distribution 
functions of the inner ES and of the additional relaxation channel remain constant, the chemical potential of the 
equilibrated outer ES is given by /i^ = -|-T(/X2 — A^i), in agreement with particle conservation within each ES and 
thereby current conservation. Finally, note that because the velocity is independent of energy, the relaxation process 
is independent of the chemical potentials of the inner ES and of the additional relaxation channel. 
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T = 0.5, T = 30 mK, Ki^t = 0.001, AE ^ 14.3 /ieV, AV = 36 /iV 
0.12n ' ' ' ' ' n 




# of iterations 



FIG. 5: (color online). Excess temperature as a function of the iteration number. Inset (a) shows the distribution functions 
after 1 iteration (i.e. the second order perturbation theory result) and inset (b) shows the distribution functions after 54 
iterations where the edge states and the additional relaxation channel have relaxed into hot Fermi functions with the same 
temperature. The dashed vertical lines in inset (b) give the asymptotic values of the chemical potential of the two edge states 
(/i°° = fj,i and fi^). 



FITTING THE THEORY TO THE EXPERIMENT 

A direct comparison of our theory with the experiment Q is difficult because the absolute values of the interaction 
matrix elements among the ES's and between the inner ES and the bulk are unknown. Furthermore, we lack the 
knowledge of the proportionality constant a between the distance and the number of iterations. Nevertheless, using 
the matrix elements and the proportionality factor a as free parameters, we could fit all of the data for different 
voltage biases and lengths simultaneously, i.e. using the same parameters to fit different values of fi2 — Mi ^-nd L. We 
find the best fit is obtained when the interaction strengths between the two outer ES's is about three times weaker 
than the interaction strength between the inner ES and the bulk excitations. The optimal value for the allowed energy 
transfer is found to be AE « 14.3 /ieV. Following the literature and taking the ES velocities to be between 10'* 
and 10^ m/s this leads to a lower bound of 0.5 /im for the momentum breaking length £p. This lower bound is smaller 
than the shortest measured propagation length of 0.8 /im and roughly one hundred times larger than the estimated 
magnetic length. 

The length dependence of the excess temperature for the optimal parameters and for different values of AV = 
— Ml, is shown in Fig. [6l where it is compared with the experiment. The thin dashed (red) curves give a least 
square exponential fit to the data. To obtain the excess temperature at an arbitrary distance x, we linearly interpolate 
the excess temperature from the numerical values obtained from the iteration procedure at distances L\_x/L\ and 
L{lx/L\ +1). This is how we obtain Fig. (3) of the main text. 

Finally, in Fig. [7] we compare the cases with and without the extra relaxation channel, using the same optimal 
parameters as before. At small voltage biases, the data agrees well with the case where relaxation occurs only 
between the inner and the outer ES's, while for larger voltage biases it agrees at long distances only with the case 
where all three relaxation channels are effective. This clearly indicates that the energy is split among three and not 
two systems in the fully equilibrated limit. 

I 
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FIG. 6; (color online). Length dependence of the 
excess temperature of the outer edge state. Com- 
parison between theory and experiment. Symbols 
with errorbars show the measurement results of 0] . 
The thin dashed (red) curves show a least square 
exponential fit to the data. 
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FIG. 7: (color online). Comparison between the 
models with (full thick curves) and without (dashed 
thick curves) additional relaxation channel. For 
small voltage biases fi2 — fii the data of 01 (black 
symbols with errorbars) is in good agreement with 
the model including relaxation among the inner 
and outer edge states only, while for larger volt- 
ages and long distances the agreement is markedly 
better with the model including an extra relaxation 
channel. 



